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We investigate hitherto unexplored regimes of probe scat- 
tering by atoms trapped in optical lattices: weak scattering by 
effectively random atomic density distributions and multiple 
scattering by arbitrary atomic distributions. Both regimes are 
predicted to exhibit a universal semicircular scattering line- 
shape for large density fluctuations, which depend on temper- 
ature and quantum statistics. 

PACS numbers: 32.80. Pj, 05.40.-a, 42.50.Lc, 71.23.-k 



Recent advances in trapping and manipulation of cold 
atoms interacting with external fields have primarily 
been implemented thus far either in single- atom systems, 
such as sparse (low-density) optical lattices 0], or in 
Bose-Einstein condensates with macroscopic numbers of 
atoms 0. 

Between these two limiting regimes lies the scarcely in- 
vestigated domain of processes involving a finite number 
of interacting atoms. Both the potential interest and the 
difficulties involved in studying 

such processes are evident in the example of optical 
lattices with appreciable atomic filling factors ||^, in 

which the transition between "insulating" (localized) 
and "metallic" (superfluid) phases have been studied in 
the framework of the Bose-Hubbard model Q]. As con- 
firmed by the above study, the mean-field approximation 
is inadequate for small numbers of interacting atoms in 
lattices due to the presence of large quantum fluctua- 
tions. Still more difficult is the analysis of fluctuations 
in systems of cold atoms coupled by long-range (1/r or 
l/r"^) field- induced interactions This leads to the 
intruiging question: Is there a way to circumvent the 
formidable task of treating the full dynamics of such 
systems and still infer their important characteristics, 
e.g., their dependence on temperature, quantum statis- 
tics (Bose or Fermi), number of atoms and lattice pa- 
rameters? And: Are there universal measurable features 
which can be a "signature" of the statistical ensemble 
(distribution function) of such systems? 

Here we consider the possibility of inferring such sta- 
tistical characteristics from the spectral features of probe 
photons or particles that are scattered by the density 
fluctuations of trapped atoms, notably in optical lat- 
tices, in two hitherto unexplored scenarios: (a) The probe 
is weakly (perturbatively) scattered by the local atomic 
density corresponding to the random occupancy of dif- 



ferent lattice sites (Fig. |l| - inset a), (b) The probe is 
multiply scattered by an arbitrary (possibly unknown a 
priori) multi-atom distribution in the lattice (Fig. |^ - in- 
set b). 

At the heart of our analysis is the idea that the Green 
function of the scattered photon or particle, which em- 
bodies the relevant spectral information, can be qualita- 
tively estimated without resorting to cumbersome pertur- 
bative calculations of the probe- multiatom interaction by 
replacing this interaction Hamiltonian by an equivalent 
random matrix. The random matrix approach, which 
has been successfully applied to various disordered sys- 
tems , allows the evaluation of probe spectra 

to all orders of scattering, expressing them by means of 
only the first two moments (the mean and variance) of the 
random interaction, averaged over the statistical ensem- 
ble of the multiatom system. The highlight of our anal- 
ysis, based on this random matrix approach, is the pre- 
diction of a semicircular spectral lineshape of the probe 
scattering in the large-fluctuation limit of trapped atomic 
ensembles. Thus far, the only known case of quasi- 
semicircular lineshapes in optical scattering has been pre- 
dicted 01 and experimentally verified |8| in dielectric mi- 
crospheres with randomly distributed internal scatterers. 

The Green function of the probe (P) at energy e = hw 
is given by 



Gp(e) = Trc 



1 



Hp-V 



PS 



(1) 



where Hp, V and ps are, respectively, 

the unperturbed probe Hamiltonian, the probe-system 
interaction Hamiltonian and the density operator for the 
ensemble of the multiatom system (S). We shall assume 
that the following conditions hold, (i) There is no appre- 
ciable back-effect of the probe on the multiatom system 
(otherwise it is no longer a probe), (ii) The state of the 
multiatom system does not change during the interaction 
time with the probe, i.e., the multiatom system remains 
"frozen", as is applicable for optical or atomic probing. 
This situation then cannot be described as Markovian re- 
laxation (exponential decay) of the probe state into the 
multiatom reservoir, since the correlation time of this 
reservoir is now much longer than that of the probe, in 
contrast with the basic assumption 

of relaxation, (iii) The probe spectrum is broadband, 
i.e., it encompasses many of its eigenstates. 
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For an ensemble "frozen" during the interaction time, 
the tracing in (|l|) implies statistical averaging over re- 
peated realizations of the multiatom system, every time 
the probe scattering is recorded, or taking the expec- 
tation value with respect to the quantum state of the 
system. For simplicity, let us explicitly consider elastic 
scattering (the extension 

to inelastic scattering is straightforward), for which 



k 



or 



(2a) 
(2b) 



Here fj: is the scattering amplitude for momentum ex- 
change hk between the probe and the system and the 
A;-mode Fourier components of the probe (system) den- 
sity operators Pp-j: (Pgj:) are defined in terms of their 
respective creation and annihilation operators Ppj: — 

E,-ata-^£, pgj: = E?-ctf'g-+fc- Equations (g) and @) 
stand, respectively, for bilinear and linear probe-system 
coupHng. For optical probes ( pa|) and (2b) correspond to 
Raman and single-photon scattering, respectively. For 
atom or neutron probes the coupling (|2^) is appropriate. 

The Green function (|l|) is obtainable, to all orders in 
V 0, by solving the set of equations for its diagonal 
elements 



k' 



where ej: are the probe energy eigenvalues in the absence 
of potential fluctuations and pointed brackets denote the 
expectation value. The spectral information contained in 
these Gj:^ is given by the density of states (DOS) of the 
probe 5(e) = ~^lmJ2kGj:^{e). 

In order to extract information on the system we shall 
make two simplifying assumptions regarding the probe 
and the coupling potential (^): (i) fj: is fiat in k (the 
coupling is strongly localized in space) within a band 
exceeding the relevant band of the system, so that w /; 
(ii) the statistical distribution of the probe is also flat in 
k and its second moment in (F^) is replacable by the 
square of its mean flux (or density) rip for the bilinear 
coupling ( pa| ) or by its mean flux (density) hp for the 
linear coupling (pb|). Under these assumptions we can 
rewrite the squared coupling potential in (||) as 



(4) 



Here (V) is the mean coupling potential and J-p ~ 
l/pnp or Tp ^ Ifl'^np in the case of ( pa| ) and (pb|), 
respectively. The quantity of interest for the system 
is the Fourier-transformed density-density correlation of 
the atomic system 



Its diagonal element 5^:^ is the static structure factor S^, 
which is the Fourier transform of the van-Hove correla- 
tion function {pg{r,t — Q)ps{r ,t = 0)) for the spatial 
density fluctuations of the "frozen" atomic ensemble. 

The difficulty of having to evaluate or measure the ma- 
trix elements Sj:-^i is avoided for a spatially random den- 
sity distribution of the atomic system, due to random 
site occupancy (Fig. |l|, inset a) and short-range interac- 
tion with the probe (e.g., a neutron or thermal atom). 
The elements S^j:i in (||) and (^ can then be replaced by 
the average of the structure factor over all relevant k: 



^kk' 



S = 



dk Sr 

k 



(6) 



where the right-hand side of S denotes the local atomic 
density or number variance averaged over the ensemble. 

The implications of evaluating the probe DOS g{e) 
using (^)-(^ will be examined for random fluctuations 
about a mean scattering potential {Vs{x)) (correspond- 
ing to the mean atomic density distribution) that is ID- 
periodic. The " unperturbed" probe dispersion associated 
with {Vs{x)) is eg = — 2J cos{kxd) + A, J being the hop- 
ping frequency, d the lattice period and A the band en- 
ergy offset. This gives rise to the following expression for 
the Green function (||) 



(3) G(e) 



k' 



(7) 



(5) 



Here {W^) = J'pS, A(e) = (W^) / y/ {e - A f - AJ^ for 
|e-yl| > 2J, A(e) = {W^)/y/A.P - [e - A)^ for |e-yl| < 
2 J and both zero otherwise. Figure [I] shows how the 
probe DOS g{e) changes from that of a periodic band 
structure corresponding to the mean potential {Vs{x)) 
to a semicircular shape as the amount of fluctuations 
measured by (W^) increases. 

In the multiple-scattering scenario, 

which pertains to resonantly scattered atomic probes 
or to intracavity optical probes (Fig. |], inset b), semicir- 
cular lineshapes are 

obtained even when the S^^i cannot be claimed to be- 
long to a random distribution (Fig. H, inset). 

In the case of strongly-interacting atoms within a lat- 
tice site or longe-range intersite density correlations |^ 
the distribution may be quite intricate, corresponding to 
sharp peaks of S^jti . Nevertheless, the universal spectral 
trends of Fig. |^ can be shown 

to hold in this scenario, provided (V^^}^^^5o(e) ^ 1, 
.go(e) denoting the "unperturbed" probe DOS. This con- 
dition allows us to estimate G^^ in (|^) to all orders in V, 
upon replacing the state of the atomic system by a gaus- 
sian random ensemble ^,0. The result is the following 
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universal formula for the renormalized probe energy e 
at a given input energy e 



(14^2 )Tr J 



(8) 



The use of (g) leads to a semicircular lineshape similar 
to the one in Fig. ^ as if the potential were random. 

In order to illustrate the role of temperature, quantum 
statistics and the mean lattice potential in producing the 
semicircular lineshape, we proceed to evaluate {W'^) for 
several simple models: 

1. The isolated-site limit : The tightly-bound Bose or 
Fermi distributions in a lattice can be estimated by tak- 
ing the potential of each site to be that of a harmonic 
well of depth Vq. The isolated-site approximation holds 
for atoms in the lowest vibrational band, when the cou- 
pling energy is much smaller than the excitation energy 



to the next band §, ^/W^ htu^ ^ y'^, A 
being the wavelength of the laser light. 
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FIG. 1. Density of states g(e) of a probe scattered by 
bosonic atoms in a ID optical lattice. Solid, dotted and 
dashed curves stand for A(e) and the thick curve stands for 
A(e) (dispersion), see text. AH curves are numerically com- 
puted from G'(e) and correspond to average random couplings 
(Vl^^) = 0.4, 2 and 10 respectively. The hopping frequency 
J = 1, and for all curves J deg{e) = 1. Inset a: A probe 
weakly scattered by a randomly occupied lattice. Inset b: A 
probe multiply scattered by a regular atomic distribution. 

In the absence of additional external perturbations, the 
coupling (W^) arises because of temperature-dependent 
fluctuations in the site-occupancy of the optical lattice, 
which has an approximately gaussian distribution 
The resulting random coupling energies averaged over 
all states yield (W^) - J^p((ng,ng,) - {nj:) {n^')) = 
JPp {{rii.nj) - {n,){nj))^. « J^p - (nsf). Here 

Tij^ = ctc^, i and j label atomic sites, and (ns) is the 
average number of atoms per site. The last step applies 
whenever rii « Uj and the density fluctuations are 
approximately site- independent. We have verified this by 
numerical simulation, considering 2 to 4 identical atoms 
on a ID lattice with 6 sites and calculating the density 
fluctuations if the probability of an occupied site is 1/10 



of the probability of an empty site. In all cases the max- 
imum relative difference between JTpiS^^, and (W^) was 
less than 10 %. 

The kinetic contribution to (VF^) due to evaporation 
of atoms from the lattice is the dominant one at high 
temperatures, regardless of the statistics. If all the atoms 
are in the lowest energy band, we may adopt the rate 
equation used to describe the formation of electron-hole 
clusters in a plasma pOj and find 



(9) 



Here a = kBtnCp, with fcp the Boltzmann constant, m the 
mass of the atoms and Cp their specific heat, T denotes 
the temperature, /3~^ = fcpT and Vq is the optical lat- 
tice potential. The influence of evaporation becomes the 
dominant effect for T ~ 25 fxK. Around T ~ 300 /iK these 
fluctuations become comparable in size to the square of 
the optical lattice potential ((VF^) ~ - 100 (neV)^) 
and atoms then largely escape from the lattice. 

At low temperatures (well below 100 fiK) the density- 
density fluctuations depend on whether the atoms in the 
lattice are bosons or fermions. For bosonic atoms in the 
lowest vibrational state we obtain jl|] 



stat,Bosc 



^ a—1 



(10) 



Here we have approximated the motion of the atoms in 
the potential wells by a harmonic oscillation with fre- 
quency uJv jlll, z = e^''^-^/''"", d denotes the average 
lattice spacing and At = {2iih'^ / mkBTY^"^ , the thermal 
wavelength, is the length scale separating quantum sta- 
tistical behavior (for At ~ d) from classical Maxwell- 
Boltzmann behavior (for Xt d). For fermionic atoms 
in an optical lattice [|3| one starts with the analog of the 
coupling (H) for particles 

with spin, using creation and annihilation operators 
ct and cg^ and performing an additional sum over the 
spin index ct, and follows the same analysis as above. One 
then finds 



stat. Fermi 



l + z 



2d3 



a=l 



.1/2' 



(11) 



with \f the Fermi wavelength. At high temperatures 
z ^ and both (jl^) and ( pd] ) reduce to the classi- 
cal Maxwell-Boltzmann result (W^^)stat,cias = z. At low 
temperatures, fermionic fluctuations approach a constant 
value, whereas bosonic fluctuations become very large as 
T decreases below l/^K, marking the Bose-Einstein 
condensation. 

In Fig. ^ we have taken typical parameters for available 
^ to show how {W"^) evolves as a 



14 



optical lattices 

function of temperature for both bosonic and fermionic 
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Li atoms. The total density-density fluctuations con- 
sist of the sum of ^ and either ( p^ ) or (|l|), depend- 
ing on the statistics. The isolated-site condition is satis- 
fed for the entire temperature range displayed in Fig. |^. 
Since the hopping frequency J ~ Vb, the random cou- 
pling for bosonic Li atoms changes from (VF^)/J^ ^ 10 
to (1^2)/ j2 _ 100, when going from T ~ 8 /^K to 
T ~ 100 /iK. Simultaneously the DOS then evolves from 
the periodic to the semicircular shape as in Fig. |l|. The 
behavior for Na or Rb atoms is found to be similar to 
that of bosonic Li atoms , apart from their different mass 
values in and (^0|). Employing currently achievable 
Bragg scattering techniques ||lj] with a far off-resonant 
laser (1 mW/cm^ intensity, 5.2 MHz detuning) and a lat- 
tice with a sufhciently high atomic filling factor the 
scattering spectrum is expected to evolve with T, in the 
microkelvin range, from the discrete to the semicircular 
regime. 
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FIG. 2. Density-density fluctuations (H^^) (in units of 10 
(neV)'^) as a function of temperature T for bosonic and 
fermionic Li atoms in an optical lattice. Thin solid line - fluc- 
tuations due to evaporation (^) (scaled by a factor of 7.8), thin 
dashed line - statistical fluctuations ( [lo[ ) . Thick solid line - to- 
tal fluctuations (W^) for bosonic Li atoms; thick dashed line 
- their counterpart for fermionic Li atoms. Parameters used 
are: Vo = 5 neV, (ns> = 0.1, d = 0.1 ^im, Cp(Li) ~ (3.6) • 10*^ 
J kg-i K-\ \f (Li) = 6 • 10"^° m and a;4Li)~ 2 • 10^ s'^ 
[ p^ . Inset: Static structure factor vs. k for phonons (solid 
curve) and nearly-free fermions (dashed curve) in a lattice at 
finite T. 

2. The nearly- free limit: A Bose or Fermi gas 
weakly modulated by the lattice potential yields S^-^i — 

^K,=k-k' ~ I'T^Kp'^i^'^^^'- Here k is a reciprocal lattice vec- 
tor, is the corresponding Fourier harmonic of the lat- 
tice potential (normalized to 1) and S^^^'^^ is the struc- 
ture factor for momentum transfer hn in a free Bose or 
Fermi gas. For a Fermi gas S^^^'^^ = ±Q{kf — k), the 
Fourier transform of pair correlations with parallel or 
anti-parallel spins (which determines the sign): it is the 
well-known step function which vanishes for k larger than 



the Fermi wavevector kf. At finite temperatures this dis- 
tribution broadens. The replacement of the nearly-free 
fermionic S^p by the average value (||) is then justifi- 
able only in the multiple-scattering scenario, while in the 
weak-scattering scenario the lattice potential harmonics 
pick out well-defined Sj:_j:,^^ (Fig. |2[ inset - dashed 
line). 

3. The phonon regime: Excitations at frequencies below 
the chemical potential of a Bose condensate trapped in a 
lattice can produce collective phonon modes |l^] whose 
"frozen" spectrum is characterized by Sg = ~^ 

1) Eg -q-G) + {nq) '}2g5{K + q + G)], where (n^) 
is the mean number of phonons at temperature T with 
wavevector q, and G denotes the reciprocal lattice vector. 
The phonon mode spectrum includes quasi-local modes 
in the case of fluctuating atomic distributions. This nat- 
urally leads to the limit and an effectively random 
coupling (Fig. ^, inset - solid line). 

To conclude, we have identified novel regimes of probe 
scattering by atoms trapped in optical lattices in the 
random-density and multiple-scattering regimes. These 
regimes cannot be treated by the mean-field approxi- 
mation, but are characterized by a universal feature of 
large density fluctuations, namely, semicircular scatter- 
ing lineshapes. This is the atom-optical analog of the 
semicircular broadening of the DOS in disordered elec- 
tronic systems, which as far as we know has not yet 
been observed unambiguously. The observation of this 
atom-optical counterpart presents a nontrivial but feasi- 
ble challenge for experimentalists. 
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